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Abstract. The Cosmic Infrared Background (CIB) peaks in the Far-Infrared 
(FIR), and its Spectral Energy Distribution (SED) is now well constrained. 
Thanks to recent facilities and Spitzer, the populations contributing to the CIB 
are being characterized: the dominant galaxy contributions to the FIR CIB are 
Luminous IR galaxies (LIRGs) at 0.5 < z < 1.5 and, to the submm CIB, Ultra- 
LIRGs at z > 2. These populations of galaxies experience very high rates of 
evolution with redshift. Because of confusion, the CIB is (and will remain in 
some domains) partially resolved and its contributing galaxies SEDs are not well 
constrained. We discuss all these aspects and show how confusion limits Spitzer 
observations, and how to overcome it in order to study the unresolved part of 
the CIB. 


1. Introduction 

A Cosmic Background that would trace the peak of the star-formation and 
metal production in galaxy assembly has long been predicted. Its discovery 
in 1996 by Puget et al. (see the review by Hauser k. Dwek 2001), together 
with recent cosmological surveys in the infrared (IR) and submillimeter, has 
revolutionized our view on star formation at high redshift. It has become clear 
that a population of galaxies radiating most of their power in the far-IR (the so- 
called “IR galaxies”) contributes an important part of the whole galaxy build-up 
in the Universe. Since 1996, detailed (and often painful) investigations of the 
high-redshift IR galaxies have resulted in spectacular results that are reviewed 
in Lagache et al. (2005). In this paper we detail briefly the sources of the 
Cosmic IR Background (CIB, Sect. 2) and their average evolution properties 
(Sect. 3). We discnss the extragalactic confusion noise, one of the most limiting 
factor of Far-IR cosmological surveys (Sect. 4) and give first hints on how we 
can overcome it with Spitzer and other facilities (Sect. 5). 


2. Sources of the Cosmic IR Background: Redshift Distributions 

The full Cosmic Extragalactic Background Spectral Energy Distribution (SED) 
is shown in figure ^ Eor clarity, we only plotted the most recent and strongly 
constraining measurements. The figure clearly shows that the optical and the 
IR cosmic backgrounds are well separated. The first surprising result is that the 
power of the IR part (the CIB, Cosmic IR Background) is comparable to the 
power of the optical one, although we know that locally the IR output of galaxies 
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Figure 1. The extragalactic background over 3 decades of frequencies from 
the ultraviolet to millimeter wavelengths. Only strongly constraining mea¬ 
surements have been reported. We show for comparison in gray an SED of 
M82 (Chanial, 2003), a starburst galaxy at L=3 10^° Lq, normalized to the 
peak of the GIB at 140 ^m. CIB measurements are, by decreasing frequency: 
Armand et al. (1994) at 2000 A; Bernstein et al. (2002) at 300, 555 and 814 
nm using the HST and Las Campanas Observatory; Mastumoto et al. (2004) 
between 2.2 and 4 p,m using the IRTS; Gorjian et al. (2000) at 2.2 and 3.3 /im 
using DIRBE and Lick] Wright (2001) at 1.25 and 2.2 /rm using DIRBE and 
2MASS] Renault et al. (2001) upper limits at 10 and 15 ^m using the CAT in 
the 7 -rays; Elbaz et al. (1999) lower limit at 15 /rm using galaxy counts with 
ISOCAM] Papovich et al. (2004) lower limit at 24 /rm using galaxy counts 
with Spitzer/MIPS] An estimate of the GIB at 60 /im from Miville-Deschenes 
et al. (2002) using IRAS] Renault et al. (2001) at 100 /im using CAT and 
DIRBE] Lagache et al. (2000) at 140 and 240 /im using DIRBE and WHAM; 
Hauser et al. (1998) at 140 and 240 /im using DIRBE] Smail et al. (2002) 
lower limit at 850 /im using galaxy counts with SCUBA] Lagache et al. (2000) 
spectrum between 200 /im and 1.2mm using EIRAS. 


is only one third of the optical (Soifer &: Neugebauer 1991). This implies a much 
stronger evolution of the IR luminosity of IR galaxies than of the optical ones. 
A second important property to notice is the long wavelength (A > 200 /im) 
behavior: the CIB slope {By cx Gispert et al. 2000), is much shallower 
than the long wavelength spectrum of galaxies (fig. This implies that the 
millimeter CIB is not due to the millimeter emission of the galaxies making the 
bulk of the emission at the peak of the CIB (at ~ 170 /im). The implications 
in terms of energy output have been drawn by e.g. Gispert et al. (2000). The 
IR production rate per comoving unit volume 1) evolves faster between redshift 
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zero and 1 than the optical one and 2) has to stay constant at higher redshifts 
up to redshift 3 at least. 

In figure n showing the CIB and a galaxy SED, we see that contributions 
from galaxies at various redshifts are needed to fill the CIB SED shape. The 
bulk of the CIB in energy, i.e. the peak at about 150 /rm, is not yet resolved 
in individual sources, but one dominant contribution to the CIB peak can be 
inferred from the ISOCAM deep surveys. ISOCAM galaxies with a median 
redshift of ~0.7 resolve about 80% of the CIB at 15 ^m. Elbaz et ah (2002) 
separate the 15 /rm galaxies into different classes (ULIRGs, LIRGs, Starbursts, 
normal galaxies and AGNs) and extrapolate the 15 /rm fluxes to 140 /rm using 
template SEDs. A total brightness of (16 ± 5) nW m“^ sr“^ is found, that makes 
about two thirds of the CIB observed value at 140 /rm by COBE/DIRBE. Hence, 
the galaxies detected by ISOCAM are responsible for a large fraction of energy of 
the CIB. About one half of the 140 /rm CIB is due to LIRGs and about one third 
to ULIRGs. However, these ISOCAM galaxies make little contribution to the 
CIB in the millimeter and submillimeter. At those wavelengths, the CIB must be 
dominated by galaxies at rather high redshift for which the SED peak is shifted. 
This effect and the redshift contribution to the CIB are illustrated in figure 
We clearly see that the submillimeter/millimeter CIB contains information on 
the total energy output by high-redshift galaxies (z>2). This is supported by 
the observed redshift distribution of the SCUBA sources at 850 ^m that makes 
about 60% of the CIB and have a median redshift of 2.4 (Chapman et al. 2003). 


Wavelength 

20% 

50% 

80% 

15 fj,m 

0.5 

1.0 

1.3 

24 pm 

0.5 

1.3 

2.0 

70 pm 

0.5 

1.0 

1.5 

100 pm 

0.7 

1.0 

1.7 

160 pm 

0.7 

1.3 

2.0 

350 pm 

1.0 

2.0 

3.0 

850 pm 

2.0 

3.0 

4.0 

1.4 mm 

2.5 

3.5 

4.5 

2.1 mm 

2.0 

3.5 

5.0 


Table 1. Redshift at which the Cosmic IR Background is resolved at 20, 50 
or 80%. Numbers have been derived using the model of Lagache et al. (2004). 

Eifty percents of the CIB are made by galaxies at redshift below 1 at 15 
and 70 /rm, below 1.3 at 24 and 160 fim and below 2, 3 and 3.5 at 350, 840 and 
2000 frm (see also Table It is clear that, from the far-IR to the millimeter, 
the CIB at longer wavelength probes sources at higher redshifts. 

So far, the most constraining surveys (in terms of resolving the CIB) are 
done at 15, 24 and 850 /rm with ISOCAM (ISO), MIPS (Spitzer) and SCUBA 
(JCMT). The deep field observations resolve about 80, 70 and 60% of the CIB, 
respectively (Elbaz et ah, 2002; Papovich et ah, 2004; Small et ah, 2002). As 
shown above, these surveys probe the CIB in well-defined and distinct redshift 
ranges, with median redshift of 0.7 (Liang et al. 2004), ~1.1 (L. Yan, priv. 
comm.), and 2.4 (Chapman et al. 2003) at 15, 24 and 850 respectively. 
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Figure 2. Cumulative contribution by redshift (from 0.5 to 8) of galaxies to 
the Cosmic IR Background, from the model of Lagache et al. (2004). Mea¬ 
surements of the CIB are reported with the same symbols as in figure^ except 
for the FIRAS spectrum, represented in black and white dashes for visibil¬ 
ity. z<0.5: oblique shaded area; z<1.0: white area; z<1.5, 2.0, 2.5, 3.0, 4.0, 
5.0: darker gray areas; z<8: solid line and horizontal shaded region. Notice 
the non negligible contribution of 5<z<8 galaxies to the CIB, at wavelengths 
above 600 /rm. 


Such well-defined redshift range selection can be understood by looking at the 
K-corrections (defined as K{L, z) = where L^(^ 2 =o) is the rest-frame lumi¬ 

nosity). This correction is specific of the spectrum of the population considered 
at a given luminosity and redshift. Taking the standard SEDs of local normal 
starbursting galaxies (i.e. with PAHs features in the mid-IR), we observe in the 
K-correction at 15 /xm a hump at z~l associated with the coincidence of the 
6 to 9 qm aromatic features and the ISOCAM filter, and at 24 qm a hump at 
the same redshift for the 24 /xm MIPS filter associated with the 11-14 qm set 
of aromatic features. Furthermore, a second hump is expected for the MIPS 
24 /xm filter at z~2, which corresponds to the redshifted 6 to 9 /xm features 
centered on thisMIPS filter (Dole et al., 2003; Lagache et al., 2004). There is 
still a debate concerning the presence of this second hump in the 24 qm sources 
redshift distribution observations (e.g. Perez-Gonzalez et al. 2005 vs L. Yan & 
P. Choi, this volume), even if growing evidences show the existence of PAHs at 
large redshift (e.g. Elbaz et al., 2005). At longer wavelengths, in the submil¬ 
limeter and millimeter, the negative K-correction becomes very effective leading 
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to an almost constant observed flux for galaxies of the same total IR luminosity 
between redshifts 1 and 5. 


3. Infrared Galaxy Cosmic Evolution 

A remarkable property of the IR galaxies is their extremely high rates of evolu¬ 
tion with redshift, exceeding those measured for galaxies at other wavelengths, 
and comparable or larger than the evolution rates observed for quasars. Number 
counts at 15 and 24 /xm show a prominent bump peaking at about 0.4 mJy and 
0.3 mJy respectively. At the peak of the bump, the counts are more than one 
order of magnitude above the non-evolution models. At 15 gm, data require a 
combination of a (l-|-z)^ luminosity evolution and (l-|-z)^ density evolution for 
the starburst component at redshift lower than 0.9 to fit the strong evolution. 
While it has not been possible with ISOCAM to probe in detail the evolution 
of the IR luminosity function, Spitzer data at 24 ^m give for the first time tight 
constraints up to redshift 1.2 (Le Floc’h et al. 2005). A strong evolution is 
noticeable and requires a shift of the characteristic luminosity L* by a factor 
(l-l-z)^-^^'^-®. Le Floc’h et al. (2005) find that the LIRGs and ULIRGs become 
the dominant population contributing to the comoving IR energy density be¬ 
yond z~0.5-0.6 and represent 70% of the star-forming activity at z~l. Those 
findings are in good agreement with the predictions from Lagache et al. (2004). 
This model constrains in a simple way the IR luminosity function evolution 
with redshift, and fits all the existing source counts consistent with the red¬ 
shift distribution, the CIB intensity, and, for the first time, the CIB fluctuation 
observations from the mid-IR to the submillimeter range. In this model, we 
assume that IR galaxies are mostly powered by star formation and hence we use 
SEDs typical of star-forming galaxies Although some of the galaxies will have 
AGN-dominated SEDs, they are a small enough fraction that they do not affect 
the results significantly. We therefore construct ’normal’ and starburst galaxy 
template SEDs: a single form of SED is associated with each activity type and 
luminosity. We assume that the Luminosity Eunction (LE) is represented by 
these two activity types and that they evolve independently and we search for 
the form of evolution that best reproduces the existing data. An example of two 
cosmological implications of this model is that (1) the PAHs features remain 
prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to 
be dominated by 10^1 Lo to ^^3 10^^ E 0 galaxies from redshift 0.5 to 2.5. 
The excellent agreement between the model and all the available observational 
constraints makes this model a likely good representation of the average lumi¬ 
nosity function as a function of redshift and a useful tool to discuss observations 
and models (e.g. figure |31 showing the contribution by redshift to the MIPS 
24 ^m source counts). Its rather simple assumptions like the single parameter 
sequence of SEDs for starburst galaxies is certainly not accounting for some of 
the detailed recent observations (e.g. Lewis et al. 2005) but probably do not 
affect seriously the redshift evolution of the averaged properties which are what 
is modeled. 

The comoving luminosity density produced by luminous IR galaxies was 
more than 10 times larger at z~l than in the local Universe. Eor compari- 
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Figure 3. Redshift contribution to the number counts at 24 /rm from the 
model of Lagache et al. (2004, their figure 7). The dot, dash, dash-dot, dash-3 
dot, long-dash correspond to the number counts up to redshifts 0.3, 0.8, 1, 
1.3 and 2 respectively. Data are from Papovich et al. (2004). 


son, the B-band luminosity density was only 3 times larger at z=l than today. 
Such a large number density of LIRGs in the distant Universe could be caused 
by episodic and violent star formation events, superimposed to relatively small 
levels of star formation activity (Hammer et al. 2005). These events can be 
associated to major changes in the galaxy morphologies. 

The large fraction of the background resolved at 850 /xm has interesting 
consequences. It shows very directly that if the sources are at redshift larger 
than 1 (as confirmed by the redshift surveys), the IR luminosity of the sources 
dominating the background is larger than 10^^ Lq, in agreement with the pre¬ 
dictions of Lagache et al. (2004), thus making a population with a very different 
IR luminosity function than the local or even the z=I luminosity function. The 
link between this population at high z and what has been seen around redshift 
one (almost by ISO) will be done by Spitzer/MIPS observations at 24 ^um. ISO¬ 
CAM galaxies contribute to about 2/3 of the energy peak of the GIB (Elbaz et 
al. 2002). The remaining fraction is likely to be made of sources in the redshift 
range 1.5 to 2.5. The presently detected submillimeter galaxies with luminosity 
IO^^Lq have an almost constant flux between redshift 1.7 and redshift 2.5 at 
24 ^m (similar to the constant flux at 850 ;um between redshift 1 and 5). The 
MIPS 24 ^m deep surveys (e.g. Papovich et al. 2004) reach an 80% completeness 
of 80 piJy (and many sources are fainter) and thus can detect all these galaxies 
when they are starburst-dominated. Gonsidering the efficiency of MIPS to cover 
large areas of the sky with a good sensitivity, it is likely that 24 ^m surveys 
will become the most effective way to search for luminous starburst galaxies up 
to z= 2.5 and up to 3 for the most luminous ones. In the near future when a 
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Figure 4. Effects of confusion in the far-IR. Observation of a source in the 
ELAIS-Nl/FIRBACK field, in a 400x400 square arcsecond box. All plates 
have been resampled to 5 arcseconds per pixel, which oversamples the far 
IR maps but undersamples the mid-IR map. From left to right: 373 mJy 
ISOPHOT 170 /im source with about 128s of integration (FIRBACK survey, 
Dole et ah, 2001); Labels indicate the 5 arcsecond pixels; Spitzer MIPS 160 gm 
with about 16s of integration (SWIRE survey, Lonsdale et al., 2004); MIPS 
70 /im with about 80s of integration (SWIRE); MIPS 24 /xm with about 160s of 
integration (SWIRE). Notice (1) the ISO 170 gm source is marginally resolved 
with MIPS 160, and is unambiguously resolved at 70 gm and 24 gm; (2) the 
two fainter MIPS 160 gm resolved sources (bottom left) create fluctuations 
in the ISO 170 gm map that produce the confusion noise when the resolution 
is limited. 


proper census of ULIRGs up to z ~ 3 will be done, the fraction of the CIB at 
~1 mm not accounted for should give an indication of the contribution from 
sources at larger redshifts. Deep surveys around 1-2 mm are the only obvious 
tool to get most of these sources. But the limiting factors of the surveys are not 
only detector sensitivity or photon noise, there is also confusion. 


4. Confusion in the Mid- and Far-Infrared 

Even if the term “confusion” might be vague, we can nevertheless dehne the 
confusion in general as the degradation of a high spatial frequency signal due to 
a poor angular resolution (see fig. El). By degradation, we mean the measure¬ 
ment of a source is affected either by a poorer quality of photometry, or a lower 
detectability. By high frequency signal we mean small-scale observed structures 
compared to the beam, usually extragalactic point-sources. Finally, by poor an¬ 
gular resolution, we mean a large beam or a poor point spread function sampling 
compared to the studied structure, usually a set of point sources. 

Predicting or measuring confusion depends on the scientific goal of the 
measurement (e.g. Helou & Beichman, 1990; Dole et ah, 2003; Lagache et 
ah, 2003): performing an unbiased far-IR or submillimeter survey and getting 
a complete sample has different requirements than following-up in the far-IR 
an already known near-IR source to get a SED and/or a photometric redshift. 
In the former case, one has to tightly control the statistical properties of the 
whole sample; in the latter case, completeness is irrelevant, and an even low 
photometric accuracy is adequate. We thus favor the use of a qualifier, like 
“unbiased confusion” for the former case. New techniques are being developed to 



use an a priori information at shorter wavelength (e.g. 8 fiui with Spitzer/IRAC 
and 24 fim with MIPS, respectively) to infer some statistical properties (like 
source density or SED) of sources at longer wavelength (e.g. 24 or 160 fim, 
respectively), and thus to beat the unbiased confusion. 

Predicting the unbiased confusion (for instance Condon, 1974; Franceschini 
et ah, 1989; Helou & Beichman, 1990; Rieke et al. 1995; Dole et ah, 2003; 
Lagache et ah, 2003; Takeuchi & Ishii, 2004; Negrello et ah, 2004) requires the 
knowledge of at least the number counts distribution of the galaxies. In practice, 
models (validated at some point by observations) are used. Since the slope of 
the counts in a Log{N) — Log{S^) diagram is varying with the flux density S^, 
the fluctuation level of faint sources below will also vary; this fluctuation level 
gives an estimate of the unbiased confusion using a photometric criterion (Dole 
et ah, 2003; Lagache et ah, 2003). At very faint fluxes, when the background 
is almost resolved, the photometric criterion will obviously give a very small 
value for the unbiased confusion level, but the observations will be limited by 
the confusion due to the high density of faint resolved sources. Thus, another 
criterion, the source density criterion for unbiased confusion (SDC, Dole et ah, 
2003, 2004), needs to be computed and compared to the photometric criterion. 
In the IR and submillimeter range below 300 gm, the unbiased confusion is in 
general better predicted by the source density criterion for current and future 
facilities, since the angular resolution has improved (e.g. from ISO to Spitzer). 
At longer wavelengths, the photometric criterion still dominates. 

In Dole et ah (2004b), we derived the confusion limits for MIPS. Using data 
at 24, 70 and 160 gm, the source density measured by Papovich et ah (2004) 
and Dole et ah (2004a) together with the modeling of Lagache et ah (2004) has 
allowed us to derive the confusion limits for Spitzer in the mid to far IR. We 
tested the model results with a Monte Carlo simulation at 24 ^m and with a 
fluctuations analysis at all 3 wavelengths (fig. ISJ. The agreement is uniformly 
very good. We determined the 5-(T unbiased confusion limits due to extragalactic 
sources: 56/iJy, 3.2 and 40 mJy at 24, 70 and 160 gm., respectively. At 24 and 
70 gm., confusion is mostly due to the high density of resolved sources, and at 
160 gm, confusion is mainly due to faint unresolved sources. Studying the FIR 
fluctuations at this wavelength is thus a tool to constrain the nature of the faint 
galaxies, beyond the confusion limit. 

What can be done with these numbers ? Generally speaking, you should 
stop integrating down to these levels, unless you know the existence of an object 
(e.g. from shorter wavelength observations) you want extract photometry from. 
It might thus be useful to integrate deeper than these limits, but the analysis 
will require some care. For instance, a sample of detected extragalactic faint 
sources will likely be highly incomplete; this might be a problem for statistical 
studies, but this can be ignored (or has to be appropriately quantified) if one 
wants photometry of known objects. 

These numbers were derived both from the data and a model, and they ex¬ 
plicitly made the assumption of a Poisson source distribution. Source clustering 
might in principle bias these estimates. Takeuchi & Ishii (2004) derived analytic 
formulas taking into account the source clustering, and showed that this effect 


Dole, Lagache, Puget: Contributions to the CIB 


9 



0 20 40 60 80 100 120 140 

time [s] 


Figure 5. Evolution atot (resulting contribution from the confusion noise 
and instrument noise) as a function of integration time, with a fit (dash) of 
the form: + ^Ln/ = . dot-dash: constant term C. 

dot: sjAjt term. Top panel (a): 24 gm. Middle panel (b): 70 ^m. Bottom 
panel (c): 160 ^m. Notice the different scales in time (seconds) and atot 
(/iJy/arcsec^). From Dole et ah, 2004b. 


is at the ~ 10% level. 

Confusion is also created in the IR by the presence of Galactic cirrus (Low 
et ah, 1984), even if its complex structure has little power at small scales (Gau¬ 
tier et ah, 1992; Herbstmeier et ah, 1998; Miville-Deschenes et ah, 2003, Kiss 
et ah, 2001 & 2003; Ingalls et ah, 2004). Usually, IR cosmological surveys are 
performed in high Galactic latitude cirrus-free regions, with a column-density as 
low as possible. For instance, only 2% of the sky have Nhi < 1-0 x 10^^ cm“^. 
Table [U gives more details about the fraction of the sky cleaner than certain 
column-densities. A common problem arises when cosmological fields, chosen 
in the visible range, happen to fall in high Galactic cirrus contamination fields, 
avoiding any efficient IR follow-up, like the 22h VIMOS field 2217-1-0024 with a 
100 gm brightness of about 4 MJy/sr. 

A number of cryogenically-cooled space telescopes have been proposed for 
the MIR, the FIR and the submillimeter spectral ranges. Table 2 of Dole et 
ah (2004b) summarizes the main characteristics of some of these observatories. 
Herschel (Pilbratt, 2001), JWST (Gardner, 2003), SPIGA (Matsumoto, 2003) 
and SAFIR (Yorke 2002), have at least one photometric channel in common 
with MIPS. As examples, Dole et ah (2004b) focused on Herschel-PAGS at 
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N(HI)/cm-^ 

Sky Fraction 

< 1.0 X 10^^ 

2% 

< 1.25 X lO^o 

5% 

< 1.6 X lO^o 

10% 

< 2.0 X 10^0 

17% 

< 2.2 X 10^0 

20% 

< 3.0 X 10^0 

30% 

< 3.8 X 10^0 

40% 

< 5.0 X 10^0 

52% 


Table 2. Fraction of the sky with HI column-densities lower than certain 
values. These numbers are derived from the HI Burton & Hartman (1994) 
Leiden/Dwingeloo survey, assuming an all sky coverage. 


75 and 170 urn, on JWST-MIRI at 24 /im, and on SPICA and SAFIR at 24, 
70 and 160 /xm, assuming in each case that the MIPS filters will be used. In 
Table 3 of Dole et al. (2004b), we use the confusion level given by the SDC, and 
compute the fraction of the CIB potentially resolved into sources. In the MIR, a 
significant step will be made with the 4m-class space telescope: as an example, 
SPICA would potentially resolve 98% of the CIB at 24 /xm. All (> 99%) of the 
CIB would be resolved with JWST or SAFIR (although doing so with JWST 
would require extremely long integrations). In the FIR, Herschel would resolve 
a significant fraction of the CIB at 70 and 160 ^m (resp. 93 and 58%, again 
with extremely long integrations). SAFIR will ultimately nearly resolve all of it 
(> 94%). 

Predictions of the unbiased confusion, for various telescope diameters and 
wavelengths, will be soon available on our website^. 


5. Beating the Unbiased Confusion 

Observing in the FIR and submm spectral ranges is relevant to characterize the 
galaxies responsible of the CIB at redshifts z > 1.3 (Table P) but the unbiased 
confusion usually avoids to directly probe these galaxies. What can be done 
? As already mentioned, one method consists in extrapolating the FIR and 
submm properties of the galaxies from the MIR (or from the radio) spectral 
properties, based on controlled samples having FIR and submm data. These 
controlled samples are usually small (e.g. Appleton et ah, 2004). This has 
the advantage of allowing to apply this extrapolation to large samples like the 
new IRAC or MIPS catalogs from GTO or legacy surveys, but the disadvantage 
of being dependent on the extrapolation based on a small sample, and of not 
allowing any object by object study. More observations are thus required in the 
FIR and submm to extend these controlled samples; more Spitzer GO programs 
should address this need in the near future. 


^http://lully.as.arizona.edu/Model 
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Figure 6 . Example of a stacking analysis at 160 ^m. 24 gm sources in the 
CDFS are selected to have flux densities 80 < S '24 < 83/iJy (330 sources). 
Images at 160 gm are stacked at the positions of the 24 gm sources. There 
is a clear detection in the resulting stacked 160 gm image in the center, with 
the expected FWHM. Pixel scale here is 1.25 arcseconds/pixel. 


Another method consists in characterizing already known submm (or FIR) 
sources. Egami et al. (2004) studied VLA and SCUBA sources in the NIR and 
MIR, and Ivison et al. (2004) did the same of MAMBO 1.2mm sources and 
both showed that IRAC and MIPS quickly detect high redshift (0.5 < z < 3.5) 
massive galaxies, dominated by star formation. Although being successful, this 
method rely also on small samples. 

Finally, a powerful method consists in making biased observations: using 
shorter wavelength data (like IRAC 8.0 gm or MIPS 24 gm) or high angular 
resolution data (e.g. VLA 1.4 GHz) to select a physically homogeneous sample, 
one can extract photometry in the FIR or submm. If the targets fluxes are 
below the MIPS FIR or SCUBA submm detection limits, one can make use 
of a stacking analysis to increase the S/N of FIR or submm data. Frayer et 
al. (2004) used this technique to extract color properties of SCUBA galaxies 
at 70 ^m. Also, Sergeant et al. (2004) stacked 450 and 850 ^m images at 
the positions of IRAC 5.8 &: 8 gm sources, and detected a positive deviation. 
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allowing them to give a constrain on the contribution of IRAC sources to the 
submm CIB. 

We present here an example of such stacking analysis at 160 /rm for il¬ 
lustration. We first selected a sample of 24 |Um galaxies with flux densities 
80 < ^24 < 83/rJy in the CDFS, and found 330 sources (Papovich et ah, 2004). 
We then stacked the corresponding 160 /rm data. Figure IHl shows the result of a 
clear detection. This illustrates the power of stacking analysis to probe fainter 
galaxy population using an a priori information from shorter wavelengths. 

Not only we can probe individual galaxies deeper into the confusion in the 
MIR using an a priori information at shorter wavelength (e.g. NIR), but we 
can also probe galaxies statistically in the FIR. In this case, obviously a careful 
sample selection is required to interpret the resulting photometry. 


6. Conclusion 

The Cosmic Infrared Background, peaking in the FIR, is made of contributions 
by galaxies located at different redshifts. The FIR CIB has most contribution 
from 0.5 < 2 ; < 1.5 LIRGs, and the submm has most contributions from 2 ; > 2 
ULIRGs. The FIR properties of these IR bright galaxies are not known in 
detail, mainly because of the unbiased confusion in this spectral range. New 
Spitzer results are starting to refine our knowledge of these galaxies. Beating the 
unbiased confusion is one of the challenges of the analysis of today’s cosmological 
surveys, and early results show that many methods are promising and successful. 

Furthermore, the efficiency of MIPS 24 fim to cover large sky areas to a great 
depth is a unique tool to search high redshift ULIRGs responsible of the submm 
CIB. There might be another Spitzer legacy consisting to cover wider areas to 
pinpoint this population, which might be characterized using also submm data 
from Herschel (and Planck), either on an individual basis or statistically. 
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